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Abstract 

We solve the Jaynes-Cummings Hamiltonian with time-dependent cou- 
pling parameters under dipole and rotating-wave approximation for a 
three-dimensional (3D) photonic crystal (PC) single mode cavity with 
a sufficiently high quality (Q) factor. We then exploit the results to show 
how to create a maximally entangled state of two atoms, and how to imple- 
ment several quantum logic gates: a dual-rail Hadamard gate, a dual-rail 
NOT gate, and a SWAP gate. The atoms in all of these operations are 
syncronized, which is not the case in previous studies [1,2] in PCs. Our 
method has the potential for extension to A''-atom entanglement, universal 
quantum logic operations, and the implementation of other useful, cavity 
QED based quantum information processing tasks. 



1 Introduction 



The superiority of quantum computing over classical computation for problems 
with solutions based on the quantum Fourier transform, as well as search and 
(quantum) simulation problems has attracted increasing attention over the last 
decade. Despite promising developments in theory, however, progress in physical 
realization of quantum circuits, algorithms, and communication systems to date 
has been extremely challenging. Major model physical systems include photons 
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and nonlinear optical media, cavity QED devices, ion traps, and nuclear mag- 
netic resonance (NMR) with molecules, quantum dots, superconducting gates, 
and spins in semiconductors [3]. 

In the last century, control over electrical properties of materials led to the 
transistor revolution in electronics. Now, in this century, PCs inspired great 
interest for controUing the flow of light. PCs are periodic arrays of dielectric 
materials, which would resemble the atomic structure of semiconductors if the 
dielectric materials were scaled down to atomic dimensions. If a defect is intro- 
duced by disturbing the periodicity of the crystal, a localized mode is allowed in 
that region. Similar to electronic band structure in semiconductors, acceptor or 
donor types of defect states can also be produced in a PBG, by locally removing 
or adding extra dielectric material, respectively. By tuning the parameters of 
the PC cavity, the frequency or the spatial profile of this localized state can be 
tailored on demand. 

Known also as photonic band gap (PBG) materials, PCs seem to be especially 
promising as a base medium both for future laser based photonic integrated 
circuits and perhaps later on for advanced quantum network technologies as 
well [4-6]. Combining a high Q factor and an extremely small mode volume 
successfully in microcavities, PCs have already become an especially attrac- 
tive paradigm for quantum information processing experiments in cavity QED 
[1,2,7]. 

Quantum logic gates and quantum entanglement constitute two building blocks, 
among others, of more sophisticated quantum circuits and communication pro- 
tocols. The latter also allows us to test basic postulates of quantum mechanics. 

In this work we propose an alternative method for creating entanglement and 
implementing certain quantum logic gates based on single mode PC microcavi- 
ties. Our scheme requires high-quality PCs with Q-factor of around 10^. Two- 
dimensional PC slabs, which have been demonstrated recently, have a Q-factor 
of only about 10^-10^. Since today's planar PC technology does not allow suf- 
ficiently high-Q cavities for our purposes, we must consider three-dimensional 
PCs despite the relative difficulty of their fabrication. 

2 Theory 

In this paper we first explore the possibility of mutually entangling two Rb atoms 
by exploiting their interaction, mediated by a single defect mode confined to a 
three-dimensional PC cavity. We then describe various logic gates that can be 
implemented using the same interaction. We achieve this in three steps: analysis 
of a single mode cavity (i) in a generic 3D PC, (ii) in a specific 2D PC and (iii) 
ultimately in a 3D version of the 2D PC of (ii). We assume that the defect 
frequency is resonant with the Rb atoms. Because the atoms are moving, the 
time-dependent Hamiltonian (Jaynes-Cummings Model) for this interaction in 
the dipole and rotating wave approximations is 
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H{t) = ^Y,^i+ hwa^a + hY^{Gj {t)a\a + h.c.) (1) 

3 3 

where the summation is over two atoms, A and B, a; is the resonant frequency, 
Gz is the z-component of the Pauli spin operator and a± are atomic raising and 
lowering operators, a and are photon destruction and construction operators, 
respectively. The time-dependent coupling parameters can be expressed as [1] 

G,(i)=no/,(i)cos(0) (2) 

where fio is the peak atomic Rabi frequency over the defect mode and is 
the spatial profile of the defect state observed by the atom j at time t. (j is 
the angle between the atomic dipole moment vector, /x^g, for atom j and mode 
polarization at the atom location. 

We ignore the resonant dipole-dipole interaction (RDDI), because the atoms 
have their transition frequencies close to the center of a wide PBG and the dis- 
tance between them is always sufficiently large that RDDI effect is not significant 
[2,8]. 

Initiahy we prepare one of the atoms. A, in the excited state and the cavity is 
left in its vacuum state, so 



l*(0)) = |100). (3) 

The PC should be designed to ahow the atoms to go through the defect. This 
can be achieved by injecting atoms through the void regions of a PC with a 
defect state of acceptor type. Since the spontaneous emission of a photon from 
the atoms is suppressed in the periodic region of the crystal, no significant 
interaction occurs outside the cavity. When the atoms enter the cavity, the 
interaction between the atoms is enhanced by the single mode cavity. This atom- 
photon-atom interaction ahows us to design an entanglement process between 
atoms. 

Given the initial state (3), the state of the system at time t should be in the 
form 



= a(i)|100) + 6(i)|010) + 7(t)|001) (4) 

to satisfy tlio probability and energy conservation. Wc will show analytically 
that the amplitudes at time t can be expressed in terms of the coupling param- 
eters and the velocities of the atoms. 

We can write the Schrodinger equation for the time-evolution operator in the 
form [9] 

ih^U{t,to) = H{t)U{t,to) (5) 
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In the basis {|100), |010), |001)}, the matrix elements of the Hamiltonian (1) are 



Hii = Hi2 = H21 = H22 = H33 = (6) 

^13 = ^^31 = flGAit) (7) 
H23 = H32 = hGB{t) (8) 

Hamiltonian operator, H{t), at different times commute, if GB{f) is a multiple of 
GA{t) by a factor of p. This condition can be easily satisfied by the appropriate 
orientation of atomic dielectric moment of the incoming atoms with respect to 
electric field in equation (2). Then the formal solution to (5) becomes 

U(t,to) = e-^'f'a'^^"'^^'> = e-TJ, (9) 

where / is defined as the integral of the Hamiltonian operator. By expanding 
the exponential (9) we obtain 

u{t, to) = i + + l^iY^'i' + - + (;^)( y)"^" + - (10) 

Multiplying equation (10) by the initial state |100) from the right and comparing 
with (4) gives 

a{t) = 1 + Gl ^(-1)" A.(G^ + Glr-' (11) 



n=l 



b{t) = GaGb E(-1)"^(G^ + (^bT~' (12) 



2n! 



7(t) = iGA E(-l)"7^;^(^^ + ^bT-\ (13) 

n=l ^ 



where we have defined 



Gj = f Gj{T)dT. (14) 

J to 

Recognizing the Taylor series for sine and cosine allows us to rewrite equations 
(11)-(13) as 

'^(*) = 1 + + G%y/' - 1] (15) 
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Figure 1: Probability amplitudes as a function of the velocity and p (i.e., red 
positive, blue negative), (a) Surface a{V,p) and (b) b{y,p) if the initial state is 
1 100). (c) Surface a{V,p) and (d) b{V,p) if the system is initially prepared in 
the |010) state. 



m = 7S^[cos(Gi + GlY'^ - 1] (16) 

7W = - \g2^ + g|)i/2 ^^<Ga + GlY'^ (17) 

For a single initially excited atom. A, passing across the cavity we can set 
Gb = in equations (15)-(17), which gives the same result with equation (15) 
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Figure 2: Slices from each surface in Fig. 1. (a) Probability amplitudes a{V,p) — 
red — and b{V,p) — dashed blue — with V = 433m/s. Entangled state (19) is 
obtained at p = 0.414. (b) Probability amplitudes a{V,p) — dashed blue — and 
b{V,p) — red — with the same velocity, V = 433m/s. Entangled state (20) is 
observed at the same, p = 0.414, value. 



of [4] . The exact solution for the time- independent problem of N identical two 
level atoms with a resonant single mode quantized field given in ref. [10] could 
be helpful to generalize our results to A^-atom case, see for example ref. [2]. 

Absent a rigorous calculation of the defect mode in the three-dimensional PC, 
let us just assume a generic spatial profile for the mode, which oscillates and 
decays exponentially. Thus, fj{t) in (2) can be expressed as [1] 

f,{t) = e~-^ cos[j{V,t - L)] (18) 

where Vj, L, Rdef, I are the velocity of atom j, the total path length of the 
atoms, defect radius, and the lattice constant of the PC, respectively. 

We choose the velocity of atoms, Vj = V, such that 150m/s < V < 650m/s, 
a typical velocity range appropriate for both experiments and our calculations. 
This has two immediate advantages: First, we make Hamiltonians at different 
times commute, and thus greatly simplify the analytical analysis. Second, we 
syncronize the atoms, providing cyclical readout that could be synchronized 
with the cycle time of a quantum computer [24]. These conditions in previous 
studies [1,2] are not satisfied. 

Using equations (15)-(18) we can express the asymptotic probability amplitudes, 
a{t) and b{t), as functions of the V and the p. The results for a{V,p) and b{V,p), 
when the initial state is |100), are displayed in Figures la and lb, respectively. 
If the initial state is |010), corresponding probability amplitudes are illustrated 
in Figures Ic and Id. Note that the surfaces in Figures lb and Ic are the same. 
This is simply due to the symmetry in equation (16). To compute these surfaces 
we used the asymptotic (constant) values of the Gj . 
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CouplingPaiameters vs. Tune 




Figure 3: Coupling parameters in the reference frame of moving atoms with 
velocities V = 433m/s at p = 0.414. VLq = llxlO^Hz, = 2.4 x IQi-'^Hz, 
I = 1.6—, L = 101, Rdef = I [!]• Atom A experiences the coupling parameter 
shown with the red solid curve and atom B the one shown with dashed blue. 



Fig. 2a illustrates a slice from the surface in Fig. la where the velocity of the 
atoms are, V — 433m/s. Note that we obtain the maximally entangled state, 

up to an overall phase, —1, when the velocities of the atoms, V = 433m/s and 
the initial state is |100). Similarly if we keep the velocities of the atoms fixed 
but set the atom B to be excited initially (i.e., |010)), we obtain the slice in Fig. 
2b and thus the maximally entangled state, 

IvPoi) - ' S (20) 

up to the same overall phase factor as in equation (19). (In equations (19) and 
(20), and in the following, we omit the cavity state in the kets, since it factors 
out and does not contribute to logic operations with which we are concerned.) 
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Piobability vs. Time 
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Figure 4: Time evolution of the probabilities that shows entanglement when (a) 
the initial state is |100) and (b) |010). 
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Figure 5: Time evolution of the probabilities that leads to a dual-rail NOT 
(PauH ax) logic operation when (a) the initial state is |100) and (b) |010). 



The calculated coupling parameters in the reference frame of the atoms as a 
function of time, Gj{t), with the specific values of fto, L, Rdef, and / are shown 
in Fig. 3. The solid red curve and blue dashed curve correspond to atom A and 
B, respectively. The total interaction time is less than 20ns. 

Probabilities |a(t)p, and |7(i)P from equations (15)-(17) are graphed in 

Fig. 4. Fig. 4a shows the evolution of the probabilities when the initial state 
is 1 100). Note that the cavity is disentangled from the atoms and we end up 
with the final state (19). On the other hand if the initial state is |010), the 
time evolution of the probabilities is illustrated in Fig. 4b. Note that the final 
state in this setting becomes (20), since the cavity is again disentangled. From 
(19) and (20) it is clear that the quantum system we have described not only 
entangles the atoms but also operates as a dual-rail Hadamard gate [3] up to an 
overall phase. 

Using the same quantum system one can also build a dual-rail NOT gate under 
certain conditions. If we set Va — Vb = 565m/s, for example, we obtain the 
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following logical transformations, up to an unimportant global phase factor, 
which defines a dual-rail NOT gate (see fig. 5): 



|10) ^ |01) 



(21) 



|01) ^ |10) 



(22) 



Furthermore using the Hamiltonian (1) it can be shown that |00) initial state 
only gain a deterministic phase factor of —1 in interaction picture. Once the 
conditions (21) and (22) are satisfied |11) initial state is transformed into itself 
up to the same global phase with the states in (21) and (22). Thus, including 
these as possible initial states our dual-rail NOT gate also operates as a SWAP 
gate up to a relative phase of the |00) state. During our analj^sis we have also 
observed that a dual-rail Z gate is also possible for some specifications. 

3 Two-dimensional Photonic Crystal 

In the preceding analysis we have assumed the generic form (18) for the spa- 
tial profile of the defect mode in order to demonstrate that, in principle, PC 
microcavities can be used as entanglers, and more specifically, as certain logic 
gates. In the following we apply these ideas to real two- and three-dimensional 
photonic crystal microcavity designs, to show that implementations of these 
quantum devices are indeed possible in these photonic systems. As the authors 
of [2] observe, however, "a rigorous calculation of the electromagnetic field in 
the presence of a defect in a 3D photonic crystal can be a difficult task". In the 
following we address this task systematically. 

First we consider a 2D photonic crystal design with a triangular lattice of di- 
electric rods with dielectric constant of 12 (i.e., silicon) in Fig. 6. The radius 
of the rods is 0.175Z, where I is the lattice constant. The symmetry is broken 
in the center by introducing a defect with reduced rod-radius of 0.0711 to form 
the microcavity of our quantum system. We assume that the atoms A and B 
travel along the dashed green lines shown in Fig. 6a, although any two of the 
obvious paths shown by dashed lines (or any void regions with line of sight of 
the cavity) work as well. 

In order to find the probability amplitudes for the atoms as a function of time 
while travelling through the crystal, we need to calculate the coupling parame- 
ters in equation (14) and substitute them in equations (15)-(17). Once we obtain 
the probability amplitudes, we can demonstrate entanglement creation and con- 
struct logic gates as before. We also note that the quality factor of the cavity 
should be high enough that the interaction time (i.e., the logic operation time) 
is much less than the photon lifetime in the cavity. Below we show that typical 
logic operations take 50/iS for the Hadamard gate and 30/Lts for the NOT or the 
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Figure 6: (a) A single-mode microcavity in a 2D photonic crystal with a trian- 
gular lattice (see the text for details), (b) Corresponding electric field spatial 
profile for the transverse-magnetic (TM) mode allowed in the cavity. 



SWAP gate. Thus a quality factor of 10^ should be sufficient for reliable gate 
operations. The quality factor of the cavity could be increased exponentially 
with additional period of rods in Fig. 6a [11]. We observe, however, that the 
spatial profile of the mode in Fig. 6b does not change significantly (and hence 
degrade the gate) after a certain number of periods. Thus, the exact number of 
periods required for a quality factor of 10^ is not essential to demonstrate our 
main goal in this paper, namely that such logic operations and entanglement 
creation are indeed possible in these photonic crystal structures. 

In our generic profile (18), the coupHng parameter is assumed to be real. How- 
ever, for generality in real photonic crystals, we must allow it to be a complex 
parameter. Thus, the interaction part of the Hamiltonian (1) for a single-atom 
cavity interaction can be written as [12-15] 

Hi ^ h\g{r)\[a''a- +ac7+] (23) 

by incorporating complex coupling parameter g(r) into (1). 

In a photonic crystal we can express the atom-field coupling parameter [12, 16] 
at the position of atom j, as 



=.9o*(rj)cos(0) (24) 
where go and '^(jj) are defined as 



- Meg / .„(,s 



>' mode 
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Figure 7: Normalized coupling parameters (i.e., divided by 50) in the reference 
frame of moving atoms with velocities V = 374m/s at p = 0.414, where go is 
found to be 2.765MHz. Blue and red curves correspond to normalized coupling 
strengths for atoms A and B, respectively. 



*(r,) ^ E{r,)/\E{r„,)\ (26) 

Vm denotes the position in the dielectric where e(r)|i?(r)p is maximum and Em 
is defined as the dielectric constant at that point. The cavity mode volume, 
Vmode, is given as 

_ JJJe{r)\E{rrdr 

Using the block iterative plane-wave expansion method [17] we found the nor- 
malized frequency of the cavity mode shown in Fig. 6b to be 0.3733c/L By 
setting I — 2.202mm, we tune the resonant wavelength to 5.9mm. At this wave- 
length, fieg for the Rb atom is 2 x 10~^^C-m [1]. Since we observe that the energy 
density is concentrated in the center of the cavity, Sm in equations (25) and (27) 
becomes 12. On the other hand, cos(Cj) in equation (24) can be safely assumed 
constant for each atom, because our cavity mode is transverse-magnetic (TM) 
mode and its electric field polarization direction can always be assumed to have 
a constant angle with the atomic dipole moment vector, fieg-' ■ We determined 
the coupling parameters as functions of time in the reference frame of moving 
atoms, as shown in Fig. 7. The tails of the coupling parameter functions do 
not contribute significantly to results due to exponential decay away from the 
cavity. 
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ibility Amplitude vs Time 




Figure 8: Probability amplitudes for entangled atoms created by a dual-rail 
Hadamard operation in the 2D photonic crystal (see Fig. 6) when (a) atom A 
is initially in the excited state and (b) when atom B is initially in the excited 
state. a(<), b(t) and j{t) are probability amplitudes for the states |100), |010) 
and |001), respectively. 
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Figure 9: Probability amplitudes for atom A and B under the dual-rail NOT 
operation in the 2D photonic crystal when initially (a) atom A is excited only 
and (b) atom B is excited only. 



With these settings Fig. 8 demonstrates an entangler which operates as a dual- 
rail Hadamard gate. If the initial state is |10) (i.e., atom A is excited), we 
end up with state l^'io) [see equation (19) and Fig. 8a]. If atom B is initially 
excited (i.e., the initial state is |01)), however, we obtain the state |^oi), up to 
an unimportant global phase factor of —1 [see equation (20) and Fig. 8b] as the 
output of the logic gate. 

In order to get the system to act as a dual-rail NOT gate, we simply set the 
velocities of atoms Va = Vb — 490m/s. The evolution of the probability ampli- 
tudes for the atoms is shown in Fig. 9. When the excitation is initially at atom 
A (or at atom B), it is transferred to atom B (or atom A). Note that we also 
get an unimportant phase factor of —1 in the output. Furthermore, as explained 
above, the input states |00) (or |11)) are not transformed into different states, 
but only |00) gain a different phase factor of H-l. Once this phase problem with 
the |00) state is solved, the system could also be exploited as a SWAP gate. 
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4 Three-dimensional Photonic Crystal 



Although the implementation of logic gates in 2D photonic crystals looks promis- 
ing, in reality we need 3D devices. The 2D analysis is useful, however, for 
reducing the substantial computation required for analyzing more reaHstic 3D 
structures, because there is great similarity between the modes allowed in these 
2D crystals and in their carefully chosen 3D counterparts (see Fig. 10), which 
we describe next. 

To demonstrate logic operations in 3D photonic crystals, we designed the struc- 
ture [18-19] shown in Fig. 10a. The 3D photonic crystal we have chosen has 
various advantages over others [20-22]: emulation of 2D properties in 3D [19], 
polarization of the modes, and simplified design and simulation. It consists of 
alternating layers of a triangular lattice of air holes and a triangular lattice of 
dielectric rods, where the centers of the holes are stacked along the [111] direc- 
tion of the face-centered cubic (fee) lattice. The parameters for the crystal are 
given in the caption of Fig. 10a. We calculated that the structure exhibits a 3D 
band gap of over 20%, across the frequency range of 0.507c/Z - 0.623c/Z. 

As in the 2D geometrj', we introduce a defect by reducing the radius of a rod in- 
side the crystal as shown in Fig. 10b. Using the supercell method we computed 
that this cavity supports only a single mode with a frequency of 0.539c/L Set- 
ting I = 3.18mm tunes the cavity mode to the atomic transition wavelength of 
5.9mm. Thus we can design a 3D single-mode cavity for our system to operate 
as the desired quantum logic gates. 

The spatial profile for the electric field of the mode is shown in Figs. 10c and lOd. 
Note that it has similar profile with its 2D counterpart in Fig. 6b, where the 
energy density is also maximized in the center of the cavity. It is this similarity 
which simplifies our design and analysis for the 3D case. 

We can quantify [19] the TM-polarization of the mode in a plane as: 



/rf2r|E(a;;r)|2 ' ^^''> 

We compute that P is almost 0.99 in the defect plane shown in Fig. 10b for 
the spatial profile exhibited in Figs. 10c and lOd. Then it is safe to assume 
TM polarized mode in the system Hamiltonian, because this doesn't affect the 
probability amplitudes significantly in the microwave regime, for the parameters 
we have chosen. 
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Figure 10: (a) Top view of the 3D photonic crystal with fee lattiee. It consists 
of alternating layers of a triangular lattice of air holes and a triangular lattice 
of dielectric rods (for details of the structure, see refs [18-19]). The nearest- 
neighbor spacing within either a hole or rod layer is -^l, where I is the fee 
lattice constant. Hole and rod radii are 0.293/ and 0.1241, respectively. The 
thicknesses of a hole layer and a rod layer are taken to be 0.225Z and 0.354/, 
respectively. Silicon is assumed as the high-index material of dielectric constant 
12. (b) Horizontal cross-section of the crystal. Dashed lines show the obvious 
paths atoms can travel. In our simulations we assumed the path shown by the 
dashed yellow Hne. A defect is introduced by reducing the radius of the middle 
rod down to 0.050/ to hold a single mode in the cavity, (c) The real part and 
(d) imaginary part of the electric-field of the TM mode allowed in the cavity at 
a particular instant in time with the frequency of 0.539c//. The imaginary part 
is half a period later. 



The coupling parameters as function of time in the reference frame of atoms are 
shown in Fig. 11, where the velocities of atoms are set to be F = 353m/s at 
p = 0.414. 

In Fig. 12 we demonstrate the 3D version of our 2D dual-rail Hadamard gate, 
which also acts as an atomic entangler. Similarly to the 2D case, if the excitation 
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Figure 11: Coupling parameters in the reference frame of moving atoms with 
velocities V = 353m/s at p = 0.414, where go is found to be 2.899MHz. Blue 
and red curves correspond to normalized coupling strengths for atoms A and B, 
respectively. 



is on atom A, the resulting state is l^'io) [see equation (19)] as shown in Fig. 12a, 
while it is on atom B the output state is |5'oi) [see equation (20)] as displayed in 
Fig. 12b, up to an unimportant global phase of —1. Interaction between atom 
A and B is mediated by the photonic qubit when correct parameters are set. 

We set the velocity of both atoms to Va — Vb — 459m/s to obtain a dual-rail 
NOT gate in the 3D photonic crystal, up to an unimportant global phase factor, 
— 1. The probability amplitude evolution of the atoms is displayed in Fig. 13. 

The excitation of the excited atom is transferred to the ground state atom 
with the help of the photonic qubit allowed in the designed cavity. Note that 
because of the symmetry of the problem, the b{t)s in Figs. 5a, 9a and 13a are 
indeed approximately the same as the a{t)s of Figs. 5b, 9b and 13b, respectively. 
Furthermore, for the same reason as in the 2D case, analyzed above, our 3D dual- 
rail NOT gate also operates as a SWAP gate up to some deterministic phase. 
That is, 



|00) ^^ -|00) 



(29) 



|01)^|10) 
|10) ^ |01) 



(30) 
(31) 
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Figure 12: Probability amplitudes for entangled atoms under dual-rail 
Hadamard operation in the 3D photonic crystal (see Fig. 10) when (a) atom 
A is initially in the excited state and (b) when atom B is initially in the excited 
state. a{t), b{t), and j{t) are probability amplitudes for the states |100), |010) 
and |001), respectively. 



a Probability Amplitude vs. Time b Probebilitj- Amplitude vs Time 




Tme (ms] Tme (ms) 



Figure 13: Probability amplitudes for atom A and B under dual-rail NOT 
operation in the 3D photonic crystal when initially (a) atom A is excited only 
and (b) only atom B is excited. 



Ill) ^ 111) (32) 

Note also that the velocities of the atoms and go for the 2D logic gates agree more 
than 90% with those in 3D. This makes sense due to the fact that the localized 
modes allowed in 2D and 3D cavities, respectively, which we considered above, 
have more than 90% overlap in their spatial profiles [19]. Thus our results are 
also consistent with ref. [19] and this is another reason, why we first investigated 
2D case. 



5 Conclusions 

To summarize, photonic band gap materials could be especially promising as 
robust quantum circuit boards for the delicate next generation quantum com- 
puting and networking technologies. The high quality factor and extremely low 
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mode volume achieved successfully in microcavitics have already made photonic 
crystals an especially attractive paradigm for quantum information processing 
experiments in cavity QED [1,2,7]. In our paper we have extended this paradigm 
by solving analytically the Jaynes-Cummings Hamiltonian under dipole and ro- 
tating wave approximation for two syncronized two-level atoms moving in a 
photonic crystal and by applying it to produce the two maximally entangled 
states in equations (19) and (20). We have also demonstrated the design of 
quantum logic gates, including dual-rail Hadamard and NOT gates, and SWAP 
gate operations. Our proposed system is quite tolerant to calculation and/or 
fabrication errors. Furthermore, most errors left beyond the design can also 
be easily circumvented by simply adjusting the velocity or the angle between 
the atomic moment vector of atoms and the mode polarization in experiments. 
Our technique could not only be generalized to A^-atom entanglement [2] but 
also has potential for universal quantum logic gates, atom-photon entanglement 
processes, as well as the implementation of various, useful cavity QED based 
quantum information processing tasks. Wc should also mention the methodolog- 
ical result that thanks to the emulation of 2D photonic crystal cavity modes in 
3D photonic crystals [19], one can design the more sophisticated circuit first 
in 2D to reduce difficulty of the 3D computations where typically much more 
computational power is needed. 
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